% CREATEMATRIX
%
% deze functie maakt van een matrix A de juist testmatrix*.in en
% testmatrix*.out files aan voor het geval van een volle matrix

function creatematrix(A,index)
delete(['../testmatrices/testmatrix' num2str(index) '.*'])

fileID = fopen(['../testmatrices/testmatrix' num2str(index) '.in'],'a');
fprintf(fileID,'%1i \n',size(A,1));
number = 30;
for i=1:size(A,2)
    for j=1:ceil(size(A,2)/number)
        fprintf(fileID,'%20.15e ',A(i,(j-1)*number+1:min(size(A,2),j*number)));
        fprintf(fileID,'\n');
    end
end
fclose(fileID)
% out
fileID = fopen(['../testmatrices/testmatrix' num2str(index) '.out'],'a');
B = expm(A);
fprintf(fileID,'%1i \n',size(A,1));
for i=1:size(A,2)
fprintf(fileID,'%20.15e ',B(i,:));
fprintf(fileID,'\n');
end
fclose(fileID)